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o: 

We present particular solutions for the following important nonlinear second order differ- 
ential equations: modified Emden, generalized Lienard, convective Fisher, and generalized 
, Burgers-Huxley. For the latter two equations these solutions are obtained in the travel- 

ling frame. All these particular solutions are the result of extending a simple and efficient 
factorization method that we developed in Phys. Rev. E 71 (2005) 046607. 

t-h ; 

SI. Introduction 

> ■ 

The purpose of this paper is to obtain, through the factorization technique, particular 
solutions of the following type of differential equations 

u + 9 {u)u + F(u) = , (1-1) 

in : 

where the dot means the derivative D = -Sp, and g{u) and F{u) could in principle be 

arbitrary functions of u. This is a generalization of what we did in a recent paper for 

Qh' the simpler equations with g(u) = 7, where 7 is a constant parameter. ^ Factorizing 

^ r-i Eq. QT1J) means to write it in the form 

-t— > ■ 
Ctf ■ 

E: 

J> \ Performing the product of differential operators leads to the equation 

•i-H , 

' u f^uu — 6\u — 62U + (frifau = , (1'3) 

Oj ' du 



[D-<p 2 (u)][D-Mu)]u = . (1-2) 



for which one very effective way of grouping the terms is 



i) 



H +<h + -p-u J ii + fafcu = 0. (1-4) 
du 



Identifying Eqs. ()1T|) and (|l-4j) leads to the conditions 

g(u) = -U 1 + 4> 2 + ^j±u\ (1-5) 

F(u) = fafau . (1-6) 

If F(u) is a polynomial function, then g{u) will have the same order as the bigger of 
the factorizing functions 4>\ (u) and 4>2 (u) , and will also be a function of the constant 
parameters that enter in the expression of F(u). 
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In this research, we extend the method to the following cases: the modified Em- 
den equation, the generalized Lienard equation, the convective Fisher equation, and 
the generalized Burgers-Huxley equation. All of them have significant applications 
in nonlinear physics and it is quite useful to know their explicit particular solutions. 
The present work is a detailed contribution to this issue. 



§2. Modified Emden equation 



We start with the modified Emden equation with cubic nonlinearity that has been 
most recently discussed by Chandrasekhar et al 2 ^ 

u + auii + f)u 3 = . (2-1) 

1) 4>i(u) = a\\[]3u , 4>2(u) = di 1 \/ r j3u, (o&i ^ is an arbitrary constant). 
Then Eq. (|1-5|) leads to the following form of the function g(u) 

gi (u) = -^(^l±l)n. (2-2) 



GSl 

Thus we can identify a = —yffi j ov ai ± = — where we use a\ as 

a fitting parameter providing that a\ < for a > 0. Eq. (|2Tj) is now rewritten as 



u 



(2ai + a^ 1 ) uii + /3u 3 = (d - a^yffiv) (d - a x y^u) u = . (2-3) 



Therefore, the compatible first order differential equation is u — ai^fpu 2 = 0, whose 
integration gives the particular solution of Eq. (|2-3|) 

1 4 . . 

u\ = or u\ = . , (2-4) 

ai ^/p(T - T ) (a ± y/cfl-80)(T - r ) 

where To is an integration constant. 

2) <j>x{u) = ai\ff3u 2 , 4>2 (u) = a^^ffi. Then, one gets 

g 2 (u) = -v^{ai 1 +3a lU 2 ) . (2-5) 

Therefore, g 2 is quadratic being higher in order than the linear g of the modified 
Emden equation. We thus get the particular case GE = 3/3, A = of the Duffing-van 
der Pol equation (see case 3 of the next section) 



u 



v//3 (c^ 1 + 3aiu 2 ) u + (3u 3 = (d - a^v^) ( D ~ ai^pu 2 ^ u = , (2-6) 



which leads to the compatible first order differential equation u — a± ^J3 U 3 = with 
the solution 

U2 = MahF • (2 ' 7) 
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§3. Generalized Lienard equation 

Let us consider now the following generalized Lienard equation 

u + g(u)u + F 3 = , 



(3-1) 



where F 3 (u) = Au + Bu 2 + Cv?. We introduce the notation A = V B 2 — AAC, and 
assume that A 2 > holds. Then: 



1) <pi(u) = ai ^ (- B + z ^) -|- Cuj , 4>2{u) = a 1 1 (J^xj^- + uj', g(u) takes the form 



9l[U) 



(B + A) (B-A) , , n 
; ai + V — ; aT 1 + (2Cai + a^ 1 ) 



2 ' 2C 
For <?(u) = gi(u), we can factorize Eq. (|3Tj) in the form 



D - a~ 



if ( B-A) 
2C 



+ u 



D-a l [^A + Cu 



u = . 



(3-2) 



(3-3) 



Thus, from the compatible first order differential equation ii — a\ ( { -^^ + Cu)u = 0, 
the following solution is obtained 

-i 



Ui 



(B + A) 



exp 



ai 



(B + A) 



(T ~ T 0) 



C 



2) (j>\(u) = a±(A + Bu + Cv?) , 02 (u) = a 1 1 ; g(u) is of the form 

g 2 (u) = - [(aiA + a^ 1 ) + 2a\Bu + SaiCu 2 ] . 
Thus, the factorized form of the Lienard equation will be 

F 3 (u) 



[D - af) 



D - ax- 



il 



u = 



(3-4) 



(3-5) 



(3-6) 



and therefore we have to solve the equation ii — a\F 3 (u) = 0, whose solution can be 
found graphically from 



ai(r - tq) = In 



ii 



In 



2Cu + B - A\^% 



{^2Cu + B + A 



(3-7) 



3) The case B = and C = 1: Duffing-van der Pol equation 

The B = 0, C = 1 reduction of terms in Eq. (|3-lj) allows an analytic calculation 
of particular solutions for the so-called autonomous Duffing-van der Pol oscillator 
equation 3 ) 

il+ (G + Eu 2 )u + Au + u 3 = , (3-8) 

where G and E are arbitrary constant parameters. Since we want to compare our 
solutions with those of Chandrasekar et al, 3 ^ we use the second Lienard pair of 
factorizing functions 4>\(u) = a\(A + u 2 ) and faiu) = .Then 



92(u) = - (Aa\ + a x 1 + 3aiu 2 ) 



(3-9) 
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Eq. (|3-8[) is now rewritten 

u - (a x A + al 1 + 3aiu 2 ) u + Au + u 3 = [D - a^ 1 ] [D -a x (A + u 2 )] u = . (3-10) 

Therefore, the compatible first order equation it— a\(A+u 2 )u = leads by integration 
to the particular solution of Eq. (|3- 1(J|) 



l li 

Aexp[2 ai A(T - t )] \ 1/2 _ / Aexp[-^AE(r - tq)] N ' 



1 - exp[2aiA(r -t )]J I 1 - exp[-§AE(r - r )] 



(3-11) 



where the last expression is obtained from the comparison of Eqs. (|3-8f) and (|3-10|) 
that gives a\ = — # and G = A ^~ 9 ■ 

This is a more general result for the particular solution than that obtained 
through other means by Chandrasekar et a/ 3 ) that corresponds to E = /3 and A = -J?. 

§4. Convective Fisher equation 

Schonborn ei a/ 4 ) discussed the following convective Fisher equation 
du 1 d 2 u du 

where the transformation to the travelling variable r = x — vt was performed in 
the latter form. The positive parameter ju serves to tune the relative strength of 
convection. 

1) = V2at(l-u), 4> 2 {u) = V2a^ . Then g(u) = -y/2 ([ai + a^ 1 ] ~ 2a i u ) ■ 

Therefore, for this g(u), we can rewrite the ordinary differential form in Eq. (|4-1|) as 

u + 2 ^- J=(oi + a^ 1 ) + V^aiuJ u + 2u(l - u) = . (4-2) 

If we set the fitting parameter a\ = then we obtain v = ^ + ^ l . Eq. (|4-2|) is 

factorized in the following form 

D-V2a^ 1 ]\D-V2a 1 {l-u)]u = 0, (4-3) 

that provides the compatible first order equation u+uu(l — u) = 0, whose integration 
gives 

ui = (liexp^T-ro)])- 1 • (4-4) 



2) Since we are in the case of a quadratic polynomial, a second factorization means 

exchanging (f>i(u) and 4>2(u) between themselves. This leads to a convective Fisher 
equation with compatibility equation u — \' r 2a^ 1 u = 0, where now a\ = — \/2/i, 
having exponential solutions of the type 

u 2 = ±exp[-ji -1 (r - r )] . (4-5) 
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§5. Generalized Burgers-Huxley equation 



In this section we obtain particular solutions for the generalized Burgers-Huxley 
equation discussed by Wang et al 5 ^ 



; du d 2 i 



+ au u —- — = (3u(l - u 6 ){u 6 - 7) , 



du 

dt dx dx 2 

or in the variable r = x — ut 

il+(u- au s )u + (3u{\ - u 5 )(u 6 - 7) = . 
1) 4>i(u) = v73ai(l — u s ) , 4>2{u) = \fj3ai 1 (u s — 7) . Then, one gets 



91 W 



\/@ (o^r 1 - ai + Ml + 5) - a^]u 5 



(5-1) 
(5-2) 

(5-3) 



and the following identifications of the constant parameters v = —\J~j3(ai — r ya 1 1 ). 
a = — \fj3 (ai(l + 5) — ajf 1 ) . Writing Eq. (|5-2|) in factorized form 



D 



/3a7 1 (n 5 -7) 



f3 ai {l-u a ) 



u = 



the solution 



(l ± exp[-ai y/p5(r - r )] 



-1/6 



(5-4) 



(5-5) 



of the compatible first order equation u — \fj3aiu{\ — u 5 ) = is also a particular 
kink solution of Eq. (|5-2|) . It is easy to solve the second identification equation for 
a>i = o,x(a,/3, 5) leading to 



-a ± y/a 2 + 4/3(1 + 5) 



Then Eq. (|5-5|) becomes a function u = u(r; a, (3, 5), and v = v(ct, (3, 7, 8). 



(5-6) 



2) </>i(n) = \fj5e\{u & — 7), (foC") = \//^ e i 1 (1 — w* 5 )- This pair of factorizing functions 
lead to 



92 ( 



(«) = ( 7 ei - er 1 + - ei(l + 5)}u l 



(5-7) 



and the ^ and a identifications: z/ = y/ft (ei 7 — e x x ) , a = y 7 /? 1 — ei (1 + 5)) . 
Eq. (|5-2f) is then factorized in the different form 



D - V^e^\l - u S )] \D - x/Pe^u 5 - 7) 



u = . 



(5-8) 



The corresponding compatible first order equation is now u — x/Pemiu 6 - 7) = 0, 
and its integration gives a different particular solution of Eq. (|5-2|) with respect to 
that obtained for the first choice of factorizing brackets: 



u 2 



1 



1 ± exp[ei v / ^7 ( 5( r - T o)] 



1/6 



(5-9) 
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U2 is different of ui because the parameter a has changed for the second factorization. 
Solving the a identification for e\ = e\{ot,(3,8) allows to express the solution given 
by Eq. (|5-9|) in terms of the parameters of the equation, u = u(r; a, f3, 7, 5), and also 
one gets v = v(a, /3, 7, S). If we set 5 = 1 in Eq. (|5-9|l . then from a = y^e^ 1 — 2ei) 

one can get e\ ± = a=l= "^g +8 — that can be used to obtain v± = v{a : j3, 7). The 
solutions given by Eqs. (|5-5|) and 1)5 -BJI and in (|5-9j) have been obtained previously 
by Wang el al 5 > by a different procedure. 

§6. Conclusion 

In this paper, the efficient factorization scheme that we proposed in a previous 
study 1 ) has been applied to more complicated second order nonlinear differential 
equations. Exact particular solutions have been obtained for a number of impor- 
tant nonlinear differential equations with applications in physics and biology: the 
modified Emden equation, the generalized Lienard equation, the Duffing-van der 
Pol equation, the convective Fisher equation, and the generalized Burgers-Huxley 
equation. 
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